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We discuss microscopic origin of integrability in Seiberg- Witten theory, following mostly the results of 
[T], as well as present their certain extension and consider several explicit examples. In particular, we 
_ , discuss in more detail the theory with the only switched on higher perturbation in the ultraviolet, 

^ ' where extra explicit formulas are obtained using bosonization and elliptic uniformization of the 

(~| I spectral curve. 



Supersymmctric gauge theories have become recently an area, which allows the application of the nontrivial 



^ ■ 1 Introduction 

l> 
vn 
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. methods of modern mathematical physics. In particular, one is often interested in the properties of the low- 
energy effective actions, which in the theories with extended supersymmetry can be expressed in terms of 
the holomorphic functions on moduli spaces of vacua, or prepotentials. These prepotentials obey remarkable 
' properties, which can be shortly characterized by the fact, that they are quasiclassical tau-functions [2], and 
. the low-energy effective theory ^3, can be formulated in terms of an integrable system [1] . 

Exact form of the prepotential provides comprehensive information about the effective theory at strong 
coupling, while at weak coupling the prepotential can be expanded over the contributions of the gauge theory 
■ instantons. As often happens in conventional quantum field theory, even each term in this infinite expansion, 
' containing the integration over the non-compact instanton moduli space, is ill-defined. It turns out, however, 
that there exists a preserving supersymmetry infrared regularization, which allows to perform a computation, 
reducing it to a sum over the point-like instantons, whose contributions are parameterized by random partitions 
[5]. Moreover, it turns our that the regularized volume of the four-dimensional space time can be re- interpreted 
as a coupling constant in dual topological string theory [6J, providing a new form of the gauge/string duality. 

This duality predicts a nontrivial relation between the deformed prepotentials of A/" = 2 supersymmctric 
gauge theories and the generating functions of the Gromov- Witten classes. Similar to the latter [7|, the deformed 
prepotentials can be expressed in terms of the correlation functions in the theory of two-dimensional free 
fermions. These correlation functions can be identified with the tau-functions of integrable systems, whose 
fermionic representation is essentially different from the conventional one [8J. We shall postpone the detailed 
discussion of this issue for the full deformed prepotentials and concentrate, following [T] , their main quasiclassical 
asymptotic. 

*Based on the talks at 'Geometry and Integrability in Mathematical Physics', Moscow, May 2006; 'Quarks-2006', Repino, May 
2006; Twente conference on Lie groups, December 2006 and 'Classical and Quantum Integrable Models', Dubna, January 2007. 
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2 Preliminaries 



Free fermions 



Let us introduce, first, the main definitions and notations for the two-dimensional theory of a single free complex 
fermion the action J tpdip on a cylinder. One can expand the solutions to Dirac equation in holomorphic co- 
ordinate w £ C*: 



so that the modes after quantization satisfy the (anti)commutational relations 

The fermionic Fock space is constructed with the help of the charge M vacuum state (a Dirac sea) 



(2.1) 



r>-M 



with 



iljr\M) = 0, r > -Af, tPr\M) = 0, r < -M 
and these definitions correspond to the two-point function 

V dzdw 



(o|^(z)^H|o) = 



(2.2) 

(2.3) 
(2.4) 

(2.5) 



z — w 



More conventional "Japanese" conventions (with the integer- valued fermionic operators ipi, if)* ^ i £ Z, see e.g. 
[S]) can be got from these by 



A^l/j I, ■0r ^ V'* , 1 : M = -n 



r+ 



r+-r 



(2.6) 



It is also convenient to use the basis of the so-called partition states: for each partition k = (fci > /j2 > • • • > 
fcf^ = > . . .) one introduces the state: 



r>-M 



and defines the U{1) current as: 

J =: tPiJj := ^ J„ 



, dw 



Obviously 



[J„, V'r] = -i'r+n, 



(2.7) 



(2.8) 



(2.9) 
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Recall the bosonization rules: 

^ =: e*"^ : , =: e"^"^ : , J = id<p (2.10) 

where 

0(2)0(0) ^- log z + ... (2.11) 

and a useful fact from U{N) and permutation's group theory: the Schur-Weyl correspondence, which states 
that 

(C*)®"- R^^TZk (2.12) 

k,|k|=fc 

as Sk y. U (N) representation. Now let U = diag (ui, . . . ,Ujy) be a U{N) matrix. Then one easily gets using the 
Weyl character formula, and the bosonization rules (|2.10p . that: 

det 

TrK,C/=(7V;k|:e'^"-i*(""):|0) = su{u,, . . . ,Uj^) = (2.13) 



detu^ 



gives the (ratio of the) standard Schur functions for any partition k, a very nice review of their properties can 
be found in lOj. In particular, from this formula one derives: 



, , \ ^ mi, , , \ ^ dimi?i, , 
|M)^5:^|M;k)=5:^|Af;k) 

k k 

with 

dim_Rk TT — k-j + j — i 

= T]— =11 —■ = 



(2.14) 



-TT (^k-^)! -pr h-kj+j-i _ Ili<,<j<e^ i^-i - ^-1 + i 



I] 



(2.15) 



being the Plancherel measure. It follows from the fact that for particular values ui = . . . = Uj^ — i 

sJ^,...,^) = ^ (2.16) 



Instantons and Nekrasov's computation 

In the context of A/" = 2 supersymmetric gauge theories, one usually starts with the microscopic theory, deter- 
mined by the ultraviolet prepotential J-^uy, which can be taken perturbed by arbitrary powers of the holomorphic 
operators 

^uv = U^o + h) Tr $2 + ^ = ^^0^' + (2-17) 

fe>0 

and quadratic JFuv = \tq Tr Then one integrates out the fast modes, i.e. the perturbative fluctuations with 
momenta above certain scale /i as well as the non-perturbative modes, e.g. instantons (and fluctuations around 
them) of all sizes smaller then /i^ . The resulting effective theory has a derivative expansion in the powers 
The leading terms in the expansion are all determined, thanks to the J\f = 2 supersymmetry, by the effective 
prepotential As ^ is lowered all the way down to zero, we arrive at the infrared prepotenial J-uv ^iK- 

The supersymmetry considerations suggest that the renormalization flows of T and T proceed more or less 



3 



independently from each other. Thus one can simphfy the problem by taking the limit, tq — *■ icxd, while Jxrv 
kept fixed. In this limit the path integral is dominated by the gauge instantons. The setup of [5^ allows to 
evaluate their contribution, as well as the contribution of the fluctuations around the instantons, exactly. The 
price one pays is the introduction of extra parameters into the problem, some sort of the infrared cutoff, which 
we denote by , since it appears to be a parameter of the loop expansion in dual topological string theory 6]. 

In particular, for the so-called noncommutative U{1) theory, or the theory on a single D3 brane in the 
background, which preserves only sixteen supercharges (so that the theory on the brane has only eight super- 
charges the instanton partition function Z{a, /i, t), t = (ti, t2, . . .), can be shown to be given by the sum over 
the Young diagrams, i.e. over the partitions [SlISlfTT]: 

^/ T^ir 1 chfc+i(a, k, ?i) 

(Z^ -p ^E/. (2.18) 



where mk is the Plancherel measure (|2.15p . and the Chern polynomials ch^+i (a, k, ft) can be introduced, e.g. 
via 

(e^_e-¥)£e"("+''(2 ^chKa,k,ft) (2.19) 

If the theory has the gauge group U{N), e.g. it is realized on the stack of N fractional D3 branes, the 
corresponding partition function is given by the generalization of (j2.18p : 

Z(a, t, h) ^ Z^^^\a, t, h) (m(5, k, h)) ' (-l)l'^lexp l ^ ^ /h^+^i k, ft) ^^^^^^ 

where m(a, k, ft) is the U{N) generalization of Plancherel measure [11] and ZP°''*(a, t, ft) is the perturbative 
partition function. 

Toda chain and tau-functions 



Consider, first, the well-known formula for the tau-function of Toda molecule (or the open N-Toda chain with 

Z{t;7i\a} 
J(t;n-l|a) 



co-ordinates 9ri(t, a) = log zft-ii-i \ a) 112)' given by all principal n-minors 



Z(t;n|a)= //K(a)2exp^t,i7,(a,J = A„x« (v4- £1 •^^) ^^m) 

K:ii<...<in l,ik 

of the N X N matrix, expressed as a matrix product with Aij ~ af ^, Dij — 6ij exp(zi) Jl^i k^i ~ a^j^^, 
where 

= YtlHl{^^) - Y^tisl + ...) + ) (2.22) 
with some appropriately chosen " initial phases" z^f^ . Rewriting (j2.2ip in the form 

^(t;"i^)=E|^i^„n ^. „3, n (^^-^.f (2.23) 

we see that the sum in (|2.2ip is in fact taken over the partitions (ki > k2 > . . . > kn) with the fixed length and 

/cj = + j - n - 1, j = l,...,n (2.24) 



^This theory can also be realized at a special point on the moduli space of U{N) gauge theory with 2N — 2 fundamental 
hypermultiplets. 
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For the particular solution of the Toda chain with ~ i, one gets for (|2.23p 

This is a singular or "stringy" solution, presenting a collection of particles, moving each with a constant speed, 
proportional to its number. In KP/KdV-theory the analog is u cx x/t, a linear growing potential of Kontsevich 
model [17J, which never topples. 

Comparing (|Z25ll to (|2?2T|) . one finds that 

-<'''i-- - (fefef )^ n ^ - -si.„. n 

In the limit N ^ oo, after particular choice of the Hamiltonians ()2.19p 

gKa.)l.^.^ ^''7|7'^ (2.27) 

and renormalization of the initial phase z^^"^ , passing from summation over partitions with a fixed length £1;, — n 
to a "grand-canonical" ensemble by a sort of Fourier transform, one gets Z{t;n) Z{a,t,fi) the ()2.18|) 

partition function. By (|2.14p . it becomes equivalent to the following fermionic correlator 



k pj-iT^^ot" — — — = {M\e--^e-^^>'>o*>'^''+^e—\M) 

where the mutually commuting modes of the M^-infinity generators can be defined as 



(2.28) 



^ ~ (2.29) 

The matrix element (|2.28p is a particular non-standard fermionic representation of the tau-function, where the 
Toda times are coupled to the VF-generators (|2.29p instead of the modes of the U{1) current (|2.8p . and it has 
been discussed in [71 [13] . 

If only ti ^0 the correlator in (|2.28|) gives 



Z{a = hM,ti,0,0, ...) = (Af|e-Tre**i^V^|M) = exp 



= exp(-:|l) (2.30) 



the partition function of topological string on P^. This is the only case when summing over partitions can be 
performed straightforwardly, using the Burnside theorem 



— 2^ aimn^ = — (2.31) 

k,|k|=fe k,|k|=fc 
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Baker-Akhiezer functions 



In addition to (|2.14p . one can consider 

Vi-. e^|Af + l;0) = ^Ck|Af;k) (2.32) 

k 

with (computed by the Wick theorem and using the properties of the Schur functions (|2.13p ) 

Ck = {M,k\i,.r e^|M+ 1;0> = JJ ^-h + r-^,^M ^^^^^^ 

i=l 

where the infinite product is actually finite 

Therefore, one gets for the Baker-Akhiezer functions 

(M|e~Tr^_^e*^fe>o*'=^'=+ie^|M + 1) 



(Af|e-^e* ^^->o tkW^+ie^ \M) 



Z{a,t,h) ^(-?l2)|k| 11 



i=l 



fc>0 

with a = and the shift S{r) = {Si{r), S2{r), . . .) generated by 



Z{a, t, h) 



(2.35) 



fc>0 

In the quasiclassical asymptotic — > 0, with hr ^ z, (|2.35p gives 



§(z, a, t, ;i) - exp i I ^ ifez'' - z(log z - 1) + a log z + . . . j (2.37) 

\fc>o / 

(similar formulas in the context of fivc-dimcnsional Seiberg-Witten theory [2] were considered in [IS]). In the 
same way one can define the two-point function 

(Af l|e-Tr?/,_^.!/;_^es5:fc>o*feM"=+ie^|M + 1) 

b (r) — 



{M\e-^ei ^fe>o *'=w^'=+ie^ \M) 
1 ^ t.h'^ (f.+ M^+i-fr-M^^+M ,2Mr'(.+l)/r(.+l) z{a,t ^ 2 ■ S{r), h) 
^'ko^ + ^ r(r+i)2 ■ Z(a,t,?i) 

S(z, a, t) 
- exp 



(2.38) 
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with 

S{z,a,t) = - 2z(logz- 1) + 2alogz + ... ^2.39) 

fe>0 

The asymptotics (12.391) plays an essential role in the study of the quasiclassical solution. Formula (|2.38p can 
be also interpreted as average of the r-th Fourier mode of the "symmetrically splitted" bi-fermionic operator 
■^(C) = (we"^/^) f/j (we^/'^) , introduced in [7]. The "doubling" of the fermions and their symmetric 

splitting along the w-cylinder turn into the double covering of z-plane by the quasiclassical spectral curve. 



Bosonization 



On a small phase space {tk = with fc > 1) the tau- function of model (I2.30p can be easily computed 
exploiting the fact that it can be presented as a matrix element of the evolution operator in harmonic oscillator 
for the initial and final coherent states. Indeed, after identifying Lq = |Jq + J-i Ji = — ^M^ + -p-a^a, where 
[a, a^] = h^, on the eigenstates with Jo ^ M ~ |, one gets for (|2.28p 

Zia, tM...)= exp (^-^ j (0|e ^e^" "e^ |0) = exp (^-^^ j g = 



exp 



where independent of the charge a part (generated by the world-sheet instantons in topological string model) 
is just a kernel of the evolution operator in the holomorphic representation with fixed at the boundaries ain 
and ttout. This result is certainly exact quasiclassically, here in the "stringy normalized" Planck constant {fi^ 
instead of h). 

If the second time t2 is also switched on, the partition function 

Z(M;ii,i2,0,0,...) = (A/|e-*e*i^°+'^*^'^«e^|M) (2.41) 
is no longer described in terms of a single quasiparticlc. Now one gets 

Lo ^ -f^o = ^ nal^an = a^a + 2A^ A + ... (2.42) 

n>0 

the system of coupled oscillators [q;„,q;3„] — ti^Snrn with the quadratic Hamiltonian, perturbed by special, 
preserving the energy due to [Hq, Hj] = 0, interaction 



More strictly, 



Wo^Hi = V2{{a^fA + a^A^) + ... (2.43) 
Z{M; tiM. 0, 0, . . .) = (0|e-^^ exp ^ {t"{a)Ho + t2Hi) ef^ |0) (2.44) 



with t(a) — ^ — h t"(a) = ti + t2a. Therefore the problem reduces to the computation of the matrix 



2 

elements 

(0|e-i^ exp - {t"{a)Ho + t2Hj) |0) = (S^)%fc)! ^^'^'''^'^ = (^.45) 
for the coherent states 

an\X) = XSn,i\X) (2.46) 
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Figure 1: Prepotential, as a sum of all connected tree diagrams in the bosonic cubic field theory. Each vertex 
is weighted by t2 and each pair of external legs - by e* "^"^ . The depicted diagramms correspond literally to the 
contribution of truncated "bosonic BCS model", with the dashed lines being the (AA''^ )- "propagators" . 

with A = and A' - Jjexpt"(a), (A|A') = exp (-^^). 
Quasiclassically, the matrix element (|2.45p gives 

(0|e-*exp-^(t"(a)ffo+t2i?/)eti|0) ^2:0 exp (^-^F(t", ta)^ (2.47) 

so that for the prepotential log Z = --^^ + 0{1) one gets from (piii)) . (PTT)) 

J^ = t(a) + F(t"(a),t2) (2.48) 

Its nontrivial part F(t"(a),t2) can be presented as a sum over all connected tree diagrams (see fig. [1]) in 
the bosonic cubic field theory (|2.45p . which is encoded by appearance of the Lambert function in the exact 
quasiclassical solution. An interesting issue would be to solve this theory exactly, at least quasiclassically, which 
is already not quite obvious even for the "truncated BCS model" of two coupled oscillators, corresponding to 
the first terms in (|2.42p . (|2.43p and diagrams depicted at fig.[T] The corresponding classical system 

a ^ t"a + 2t2a^A, = -t" - 2t2aA'^ 

. .2 (2-49) 
A = 2t"yl + ^2a^ ^ -2t" - t2a^ 

can be integrated in terms of elliptic functions, and possesses, in particular, a "kink" solution. We shall return 
to a detailed discussion of these issues elsewhere. 



3 Quasiclassical free energy 

The instantonic calculation in A/" = 2 extended gauge theory (|2.17p gives rise to the Seiberg-Witten prepotential 
as a critical value of the functional! 

dxf'ix) J2 dx,dx2f"{x^)f"{x2)F{x, - X2) (3.50) 

,2 £ 

extremized w.r.t. second derivative of the profile function f"{x) — with the kernel 

F(x) = ^(^log:r-0 (3.51) 

coincides with the perturbative prepotential of pure M = 2 supersymmetric Yang-Mills theory. Formula p.50p 
means that the quasiclassical free energy for the partition functions (|2.18p and (|2.20p is saturated onto a single 



^We choose here different (by a factor of 2) normahzation of the time variables t compare to [l] and correct some misprints. 
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"large" partition k* with the profile function /k* (x) — f{x), where for each partition k = fci > ^2 > • • • > > 
fc£^+i = 0, . . . the profile function is defined by 

fkix) — \x — a\ + {\x — a — h{ki — i + 1)| — |x- — a — h{ki — i)\ — \x — a — h{l — i) \ + \x — a + hi\) (3.52) 

(see [Hill] for details). In particular, one can write for (|2.19p 

chi(a,k) = -/ dx/i;;(a;)x'-^((a + ;i(fc, -z + l))'-(a + ft(fe, -i))') (3.53) 

The variational problem for the functional (j3.50p should be solved upon normalization condition for f{x) and 
the constraint 

a = i I dx xf"{x) (3.54) 



which can be in standard way taken into account by adding it with the Lagrange multiplier 

T^T + a^[a-\ J dx xf"{x)^ (3.55) 

having a sense of the k = term in the summation in formula (|3.50p . The whole setup of (|3.50p is almost 
identical to the standard quasiclassics of the matrix models, where the Coulomb gas kernel is replaces by a 
(multivalued!) Seiberg-Witten function p.5ip . 

The extremal equation for the (|3.50p gives 



J2 ikx'' - I dif"{x){x - x) (log \x-S:\-l) = a" 

k>0 



(3.56) 



on the support I where f"{x) ^ 0. Generally, for the microscopic non-abelian theory this support consists of a 
set of several (disjoint) segments along the real axis in the complex plane, where the filling fractions are fixed 
separately with the help of several Lagrange multipliers, see below. Equation (|3.56p means that 



S{z) = J2 tkz'' - J dxf"{x){z - x) (log(z - x) - 1) - 



(3.57) 



fe>0 

is an analytic multivalued function on the double-cover of the z-plane with the following properties: 

• The real part of the multivalued function (|3.57p vanishes 

S{x) = \{S{x + iQ) + S{x-iQ))^Q, xel (3.58) 

on the cut, due to (j3.56p . 

• For its imaginary part one can write 

{0 , z>x+ 
a-z + f{z), x^<z<x+ (3.59) 
2(a — z) , z < x~ 
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• We see from p.59p that even the differential dS is muhivalued. Indeed, one can easily establish for 



dz 

that 

1 



$ = ^ = t"(z) - / dxf'ix) log(z - x) (3.60) 



TT 



, z > x^ 

lm^z±iO) = ±{ l-/'(z), x^<z<x+ (3.61) 

2 , z < x^ 



However, the differential 

d$ = t'"{z)dz + dz [ (3.62) 
J z- X 

is already single-valued on the double cover of the cut z-plane with the periods / d^ ~ Airi'L, so dS is 
defined modulo Airidz, and one can make sense of the periods / dS due to f dz = 0. Therefore, the 
exponent exp ($/2) is already single- valued on the double cover and equals to unity on the cut. 

• In order to consider the asymptotic of (|3.57p in what follows we shall always choose a branch, which is 
real along the real axis, i.e. take it at real x — s- +oo. In particular, all residues below could be understood 
in this sense, as coefhcients of expansion of generally multivalued differential at x ^ +oo. 

• Taking derivatives of ()3.56|) in x-variable, or integrating by parts, one can bring it literally to the form, 
arising in the context of matrix model. However, for the purposes of Seiberg-Witten theory one needs a 
solution with different analytic properties: in matrix models the resolvent G ^ ^ does not have poles at 
the branching points where dz — Q (see e.g. |16) and references therein), which is not true for (|3.57p . 

Asymptotically from (|3.57|) one gets 

r.k+1 



S{z) = ~2z(logz-l)+y"tkz'' +logz f dxxf"(x)-a^ -2\^ f dxf"(x) 

fc>0 fc=l 

= -2z(logz - 1) + ^ t,z^ + 2alog. - 1^ - 2f: 

fc>0 k=l 

where, according to (|3.50p . (with convention that reSoo^ — 1) 

dT _ 1 
dtk ^ 2(fc + 1) 

and, due to (|3.55|) 



x 



(3.63) 



dxf"{x)x''+^, k>0 (3.64) 



oa 



(3.65) 



The coefficient at the z(logz — 1) term is fixed by normalization 

J^dxf"{x) = f{x+)~f{x-)^2 (3.66) 

where x^ (in the one-cut case) can be defined as two solutions to the equation 

/(x±) = |x±-a| (3.67) 
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Using variational equation (|3.56p . one can also write for the functional (|3.50p the double- integral representation 
(of. with [Hj) 

^=\l dxidx2f"{xi)f"{x2)F{xi-X2)+aa^ + ^0 (3.68) 

^ Jxi>X2 

expressing it in terms of the perturbative kernel p.5ip and extremal shape f{x), solving p.56p . The (time- 
dependent) constant Sq arises in (I3.68P due to constraint (|3.66|) and appears to be the constant part of the first 
primitive of the function (|3.57p . 



4 Dispersionless Toda chain 

In the case of a single cut let us present the double cover of the z-plane = [z — a;+)(z ~ x^) in the form 

z = v + k(^ + ^ (4.69) 

with x^ = u ± 2A and 

= {z- vf - (4.70) 

On the double cover (j4.69p . which is in the case of single cut just with two marked points P±, with z[P±) = oo, 
w'^^{P±) = oo, formula p.57p defines a function with a logarithmic cut and asymptotic behavior (|3.63p . odd 
under the involution w <-> — of the curve (|4.69p . In terms of the uniformizing variable w one can globally write 



(4.71) 

fe>0 



S" = -2 + A + j logw - 2A(logA - 1) - ^ + ^tkVlk{w) + 2alogw 
where 

nk[w) ^ zX- z^L, k>0 (4.72) 

are the Laurent polynomials, odd under w ^ ^. The first term in (|4.7ip comes from the Legendre transform 
of the Seiberg-Witten differential - z-^^. 

The canonical Toda chain times are defined by the coefficients at the singular terms in (|3.63p 

to = resp_^^dS ^ -Tesp_dS ^ 2a (4.73) 

and 

tk^v resp^z-'^dS ^ ~y Teap_z-'''dS, k>0 (4.74) 
k k 

From the expansion (|3.63p it also immediately follows, that 

^ = -Yeiip^z''dS ^ --resp z'^dS, k>0 (4.75) 
dtk 2 + 2 ' ^ ' 

Formulas (|4.73p . (|4.74p . (|4.75p together with p.65p identify the generating function (I3.50p with the logarithm 
of quasiclassical tau-function, being here, in the case of a single cut, a tau-function of dispersionless Toda chain 
hierarchy. 

The consistency condition for ()4.75p is ensured by the symmetricity of second derivatives 

^'•^ ^-lesp^iz'^drin) (4.76) 



dtndtk 2 
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where the time derivatives of (|3.63p 



^ dS / , d^T ^yr^ d^T 1 



, ^ n>0 / ^4^^^^ 

dtk z^P± \ dadtk otkOtn nz'^ 

\ n>0 / 

form a basis of meromorphic fimctions with poles at the points P±, with z(P±) = 00. All time-derivatives here 
are taken at constant z. 

Expansion (|4.77p of the Hamiltonian functions (|4.72p expresses the second derivatives of in terms of the 
coefficients of the equation of the curve (|4.69p . e.g. 

Qq = 21ogz - 21ogA- — - 

z^oo z z- 

2A2 2uA2 , . 

ill = z-v ^ + . . . (4.78) 

z— ►oo z Z^ 

2 . 2 4vA2 2A2(A2 + 2f2) 
^2 = z - {v^ + 2k^) ^ '- + ... 

2— too z 

Comparison of the coefficients in (|4.78p gives, in particular, 
and 

which becomes the long-wave limit of the Toda chain equations after an extra derivative with respect to a is 
taken 

= Aexp— (4 81) 

dil da da 

for the Toda co-ordinate — Substituting expansions (|4.78p into (|4.76p . one gets the expressions for the 
so called contact terms [19] in the U{1) case, which are the polynomials of a single variable v with A-dependent 
coefficients. 

One can now find the dependence of the coefficients of the curve (|4.69p on the deformation parameters t of 
the microscopic theory by requiring dS = at the ramification points, where dz = 0. This condition avoids from 
arising of extra singularities at the branch points in the variation of dS w.r.t. moduli of the curve. Equation 

'^^ k(w~-\^Q (4.82) 



d log w \ w 

gives w = ±1, where now 



dS 



E dUk 
tk 

k>0 



-I- 2a - 2w T 4AlogA = (4.83) 



d\ogw 

If t/c = for fc > 1 , solution to (|4.83p immediately gives 

v = a, A^ = e*i (4.84) 
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and the prepotential 
Adding nonvanishing t2, one finds 



^ = + l'-^^P+ i^dS) - y = \ah, + e*^ (4.85) 



2i2 ^ (4.86) 
log = ti + 2t2a - L (-4i2e*i+2*2'') 



where the Lambert function L(t) is defined by an expansion 

L(i) = f:^P^=^-^' + ^*'-^^' + --- (4-87) 

n— 1 

and satisfies to the functional equation 

L(i)e^(*) = t (4.88) 
Hence, for the prepotential with ti , t2 7^ one gets 

= ioo" + J 11(1 - t)<t»8p, (-''<iS) + jresp^ (jdS) - y = (4.89) 

fc>l 

11 1 a?' 

, ^i.-.^ 9°°^ + oi^esp^ (zdS*) - -i2resp^ {z^dS) - — 

tfc=0,K>2 Z Z 4 Z 

wherefrom the instanton expansion can be computed (which can be strictly got as an expansion in parameter 
q after ti — > + ^ logg) 

T - t(a) + e^''^-^) + tl e2t"(-) + g3t"(a) _^ e4t"(a)_^ 
1 fin „ 1 c^'^R ,, 1 (4.90) 

^IWJ^S g5t"(a) ^ i:^^10 g6t"(a) ^ ^ ^ ^ ^ ^^^^ ^ ^ 75(a)5"(a) + • • • 

O O 4: 

with 5(a) = expt"(a). Expansion (I4.90p directly corresponds to summing over connected tree diagrams in 
bosonic model, presented at fig. [TJ 

It is also easy to compute the explicit form of the extremal shape for nonvanishing ti,t2, which reads 



fix) = ^ (^arcsin + h^J - [x - vf^ 



(4.91) 



w-2A<a;<'i; + 2A 

where v — v{t) and A = A(t) are given by (|4.86p . or obey 

'y-a = 2t2A^ log A^ = ii + 2t2f (4.92) 

Formula ()4.9ip is a direct consequence of (I3.62p and directly following from ()4.7H) upon the relations (|492l) 
expression 

.-y-.z — V a — V + 2t2A-^ 
$(w; ii, i2, a) = 2t2y - 2 logw + {h + 2<2W - log A^) + 2 = 

^ ^ (4.93) 



^ 2t2A \w 2 log w 
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Note, that (|4.9ip stays that the Vershik-Kerov "arcsin law" [2U] for the Hmiting shape is deformed by the Wigner 
semicircle distribution. 

If the first three Toda times ti,t2, is are nonvanishing, instead of (|4.9ip one gets 

fix) = ^ (^arcsin (^^) + (t2 + 3hv)^4A^ ~ {x ^ v)^ + ^^^x - v)^4A^ - {x - v)^^ 



(4.94) 



i;-2A<a:<w + 2A 

where v and A are now subjected to 

V -a^ 2{t2 + 3t3v)A^ 
log A^ = ti+ 2t2v + Msiv'^ + 2K^) 

Generally we obtain for the limit shape 



(4.95) 



f'{x) = ^ l^arcsin + Y.^tkQk{xW AK-^ ~ {x - vf^ 



v-2K<x<v + 2K 



(4.96) 



where v and A obey some sort of hodograph equations v — a = P„(v,A;t), logA^ = PA(w,A;t) for some 
polynomials Py and Pa, whose expansion in Toda times t can be easily reconstructed from the presented above 
formulas of general solution. 

5 Extended non-abelian theory 

In the case of U{N) gauge theory one has to consider solution with N cuts {!;}, i — 1 . . . , N, which arises after 
adding to the functional (|3.50p N constraints with the Lagrange multipliers 



J'^J' + J2af(^a,~^J^dxxf"{x)^ (5.97) 



i.e. solution to the integral equation 

J^tkx'' - J dxf"ix)ix-x){log\x-x\-l) = a^, xeh, i = l,...,N (5.93) 



fe>0 

Now it can be expressed in terms of the Abelian integrals on the double cover 

N 



y' = Y[iz~xi)iz-x-) {5. 



i=l 

which is a hyperelliptic curve of genus g = N — I. Define, as before: 



S{z)^t'{z)- / dxf"{x){z-x)ilog{z-x)-l)-a" (5.100) 



where the integral is taken over the whole support I — U^^Ii, = Sjli ^ ^'^'^ consider its differential, 



N ^3 = 1 "'j 

or 



"■^ k>o 



(5.101) 
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satisfying 

$(x + iO) + $(a;-iO) = 0, xeh, i=l,...,N (5.102) 



on each cut, and normalized to 



^x+) - 0, 



^{xj ±iO)^<S>{x+_j^±iO)^±2TTi{N~j + l), j = 2,...,N (5.103) 
<^ix^) = ±2'KiN 



Vanishing microscopic times 

Consider, first, all tk — for fc ^ 1, and define A^^ = 6*^. Now $ = ^ is an Abelian integral on the curve 
(|5.99p with the asymptotic 

$ T2N log z±2N log A + 0{z-^) (5.104) 

whose jumps are integer- valued due to ()5.1G3|) . or <f ^ AniZ. It means that the hyperelliptic curve ()5.99p can 
be seen also as an algebraic Riemann surface for the function w = cxp (— <i>/2), satisfying quadratic equation 



N 

\ 7/1 / 

i = l 



1 \ ^ 

-) = PNiz)^Y[{z~v,) (5.105) 



since for the two branches 1/7+ = w and w_ = one immediately finds that their product w+ ■ W- and sum 
+ ui_ are polynomials of z of given powers (zero and N correspondingly) . 

Equivalently, the ends of the cuts in ()5.99|) are restricted by N constraints in such a way, that this equation 
can be rewritten as 

^PN{zf -4A^^ (5.106) 

i.e. {xf} are roots of Pn{z) =F 2A^ = 0, and 



y^K^lw--] (5.107) 



The generating differential (|5.10ip is now 



dS = -2logwdz = ~d(2zlogw) +2z— (5.108) 

w 

just the Legendre transform of the Seiberg-Witten differential dS ~ on the curve (|5.105p . (|5.106p . It 
periods 

ai^ — i z — 5.109 

2Tn Ja, w 

coincide with the Seiberg-Witten integrals and the only nontrivial residues at infinity give 

resp_^ {z~^dS) = -resp_ (z-^dS) = logA^^ 

^ (5.110) 



resp_^ (dS) = -resp_ (dS) = 2 ^ 
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The differential (|5.108p satisfies the condition 



CO ^ SP{z) , 

dab ^ — az — az — nolomorpnic (5 m) 

w y Ef=ii'j=o 

where the variation is talcen at constant co-ordinate z and constant scale factor A. Thus, the integrable system 
on "small phase space" is solved for the scale A^^ = e*i and the moduli vj , j — 1, . . . , N oi vacua of the U{N) 
gauge theory, satisfying the equation X^^i — ^ ^-^i^ ^^l^ transcendental equations for the Seiberg-Witten 
periods (|5.122p . 



Nonvanishing microscopic times 

When we switch on "adiabatically" the higher times (|2.17p with fc > 1, the number of cuts in (|5.98p remains 
intact, and the differential (|5.10ip can be still defined on hyperelliptic curve (|5.99p . However, now the role of 
bipole differential ^ of the third kind is played by 



d$ = I ^ k{k ~ l)tkz''-^ - / 
\fe>i •' 



dxf'ix) 



Z — X I 

(5.112) 

= ^ fc(fc - l)tkdnk-i - 2Ndno - 47ri ^ dioj 
k>i j=i 



where duji, i = 1, . . . , — 1 are canonical holomorphic differentials normalized to the A-cycles, surrounding first 
A^ — 1 cuts. The differentials dflk in (|5.112p are fixed by their asymptotic at 2; oo 



kz''-^dz + 0{z-^), k>0 
+ 0(z"^), k^O 



. = „ { dz ^, , „ (5.113) 



z 

and vanishing ^-periods 

df]fc = 0, fc>0, Vi=l,...,7V-l (5.114) 



The nonvanishing periods of (i$ are fixed by 

d$ - -2m [ f"{x)dx = -2m {f{x+) - fixj)) = -4TTi (5.115) 

which justifies that generating differential dS is still defined modulo Airidz. The only important difference with 
the previous case is that integrality of the periods f d^ ^ 47riZ, which was reformulated in terms of an algebraic 
equation (|5.105p for the theory on "small phase space" , remains now a transcendental equation, which cannot 
be resolved explicitly. 

Nevertheless, on the curve (j5.99l) any odd under hyperelliptic involution differential can be always presented 

as 

d$ = f(f)^ (5.116) 

y 

where s{z) is a polynomial of power N + K — 2 m case of nonvanishing microscopic times ti, . . . , t/^ up to the 
AT-th order. Its higher K coefficients are fixed by leading asymptotic (i2, . . . , ii^) and the residue at infinity 

resp±d$ = T2N (5.117) 
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and the rest N — 1 coefficients can be determined from (|5.115l) . Tlris fixes completeiy tlie difFerentiai on tfie 
cm've (|5.99[) which still remain to be dependent upon 2N (yet arbitrary) branch points {xf}. 

The generating difFerentiai dS can be now defined in terms of the Abelian integral $(z) 



dS = ^dz ^dz / d$ (5.118) 



zo 



The dependence upon 2N + 1 parameters (the positions of the branch points in (|5.99p together with zq) is 
constrained by additional to (|5.115p vanishing of the B-periods 

I d<P = 0, j = l,...,7V-l (5.119) 

Integral representation (|5.10ip suggests a natural normalization (|5.103p . i.e. 

zo^x+, <^izo) = ^x+) = (5.120) 

where is the largest among real ramification points {x^}. These conditions lead to the following form of 
expansion of <l>(z) in the vicinity of ramification points 



^z) =^ <^{xf) + (^f^z~xf + 



2s{xf) (5-121) 



where the constants ^{x^) are given by (|5.103p . 

The rest iV + 1 parameters are eaten by the periods 

a, = i dxxf'ix) = £ zd^ = i - 1^ ■ • ■ > ^ - 1 (5.122) 

together with the residues 

a = ^ J^dxxf"{x) = -iresp^ {zd<^>) (5.123) 

and the "free term" or scaling factor 

= resp+ (z-i$dz) (5.124) 

Recall once more, that an essential difference with the case of vanishing times is that for tk ^ 0, the exponent 
exp (<I>) acquires an essential singularity at the points P±, and the constraints (|5.115p . (|5.119p cannot be resolved 
algebraically. The form of the expansion (|5.12ip ensures that variation of the generating differential at constant 
z w.r.t. moduli of the curve (|5.99p 

5{dS) = 6 {<^dz) = 

-s{xf)5xf dz (5.125) 
= ^ ~ holomorphic ^ ' 



'""'^ rife \ i^, - xt)ix, -x^) Jz-x 



is indeed holomorphic. 

The Lagrange multipliers 



dai 



(5.126) 
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can be computed by a standard trick. Consider equation (|5.98p for j and fix there a;-variables to be at the 
ends of corresponding cuts. Then 



1 

2 ./a. 



af-af = Re/ dS ^ - (j) dS (5.127) 



or 

g = ^ = l,...,N-l (5.128) 

For the time-derivatives of prepotential one can write 

6 Quasiclassical hierarchy and exphcit results 

From the expansion (|5.10ip in the case of U{N) extended theory it still follows that the first derivatives of 
quasiclassical tau- function T are given by (|3.64|) and (|4.75p , while for the second derivatives one gets (|4.76p , or 

= iresp+(z™df^„) = iresp+«p^ (z(F)"z(P')"W(^, ^')) (6-130) 

where we have introduced the bi-differential W{P,P') — dpdp' \ogE{P, P'), with E{P,P') being the prime 
form, see [21] for the definitions. In the inverse co-ordinates z = z{P) and z' — z{P') near the point P+ with 
z{P~^) = oo it has expansion 

''^ = + • • • = E iITT '^"^(^') + ■ ■ ■ (6.131) 

^ ^ k>0 



The bi-differential W{P, P') can be related with the Szego kernel 



Se{P,P')S-,{P,P') = W{P,P') + doj,{P)dLOj{P')^\ogeMT) (6.132) 



which, for an even characteristics e = — e, has an explicit expression on hyperelliptic curve (|5.99p 



5.(.,z') = i^2(f)±m^ (6.133) 
2^U,{z)U,{z') z-z' 



with 



Ue{z) 



N 



z - x + 



Here x^± is partition of the ramification points of (j5.99p into two sets, corresponding to a characteristic e. 

For example, on a small phase space, when (|5.99p turns into the Seiberg-Witten curve (|5.105p . there is a 
distinguished partition e ^ E, corresponding to an even characteristics with 

, , /P(z)-2A^ 
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Substituting (|6.132p . (|6.133p into (|6.130p gives 



a, a, i'-esp.^p^ (z(P)"z(F')"5e(F,P')') - 

-iresp^ • resp^ ^ log0e(O|T) = (6.136) 

a2 -r- a2 -r- a 
S OaiOtn OajOtm Olij 



where for the "contact polynomials" one gets from (|6.133p 



= ^resp.^p, (^^— + + 2CJ^ j 



dzdz' (6.137) 



If calculated in the vicinity of the small phase space and for the particular choice of characteristic (|6.135p . 
residues (I6.137P vanish for n,m < N, and one gets exactly the conjectured in [19 formula 



N 



u„ = 2^;^ ^ -resp^ [ ] = Ij^v? ^ -(Tr (6.139) 



dt„dtm 2 dai daj dTij 

with 

. = -resp_^ z'"-^ — 

atn-i n \ y J n 

Equation (I6.138P is a particular case of the generalized dispersionless Hirota relations for the Toda lattice, 
derived in [T5j . 

Let us point out, that this derivation of the renormalization group equation (I6.138P in [T is almost identical 
to developed previously in 22J for another version of extended Seiberg-Witten theory, which can be defined by 
generating differential 

dS ^Y.TA ^Y.TkP[z)T^ (6.140) 

fc>0 fe>0 

directly on the Seiberg-Witten curve (|5.105p (whose form remained intact by higher flows in contrast to 

the quasiclassical hierarchy, determined by (|5.118p ). and all derivatives were taken at constant w. For example, 
if = 1 and only Tq, Ti do not vanish with d^li = dVli = (z — w)^, one gets 

^ ^z-v) = ?^(w+^ (6.141) 



dTi ' ' dTi\ w 

and 



= dVLi + Ti^— w + — — = 1 + -tt;— (6.142) 



aTi dT^X w) w \ A dTi 

which means, in particular, that the scale factor A ^ Ti linearly depends on the first time, in contrast to the 
exponential dependence in formula (|4.84p . Indeed, taking derivatives of (|4.7ip at constant z, instead of (|6.14ip 
one gets 

dv ^ dA f , ^\ , I, f l \ d log w 

and therefore 



= ^ (tiVli + 2alogw - 2zlogw - 2A(logA- 1) ( w - - ) ) = A 
oti oti \ \ wj J 

+ ^(ii-21ogA) \ \ + —Q^ f (ti-21ogA)Afw+---j + 2a - 2w 



(6.144) 
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i.e. formulas ((i?72)) . (|iT71) are provided directly by 

The choice of extension (|6.140p in [5^ was motivated rather by technical reasons: preserving the form 
of the Seiberg-Witten curve (|5.105p with deforming only the generating differential, moreover that the latter 
remained single-valued even in the deformed theory. We see, however, that the old choice is not consistent 
with the microscopic instanton theory (|2.17[) . (|2.18p . (|2.20p . basically since the appropriate co-ordinate for the 
quasiclassical hierarchy is z, coming from the scalar field <i> of the vector multiplet oi M = 2 supersymmetric 
gauge theory. However, the corresponding quasiclassical hierarchy is defined even more implicitly, due to the 
highly transcendental ingredient § ^ AiriZ for the second kind (not for the third kind) Abelian differential, 
and one needs to apply special efforts to extract explicit results. 



Instanton expansion in the extended theory 

The instantonic expansion J- = X]a:>o •^'^ non- Abelian theory starts with the perturbative prepotential 

N 

defined entirely in terms of the functions (|2.17p and p.5ip . ft is totally characterized by degenerate differential 

d^o = t"'iz)dz - 2^^^ ^ t"'{z)dz - 2 y (6.f 46) 

Pn{z) ^z-Vj 

(which does not depend on higher times), and the coefficients of the polynomial Pn{z) in (|5.105p . (|6.f 46p 
coincide with the perturbative values of the Seiberg-Witten periods 

Gi = —^leSy.zd^o — Vi (6.147) 

The perturbative generating differential is dSo — ^odz, with 

JV 

$o = t"(z)-2^1og(z-t;j) (6.148) 



and satisfies 



f^ = 2^, , = 1,...,A^ 

oaj z — Vj 



(6.149) 



dtk 



kz" ^dz, k > 



what gives rise to 



Equations 



N 

So{z) = t'(z) - 2^(z - V,) (log(x - Vj) - 1) (6.150) 
«f = g^ = 2^o(a,) (6.151) 



^1 

completely determine (|6.145p . since on this stage one makes no difference between Vj and a^. 
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Moreover, vanishing of the S-periods (|5.119p of the differential (|6.146p 

(i$o = (6.152) 



where = aj ± y^qSj + 0{q^) are positions of the branching points of the curve ()5.99p in the vicinity of 
perturbative rational curve, immediately gives the deviations 

pt"(ai) 

= l,...,iV (6.153) 



where the numeric coefficient is fixed from comparison with the Seiberg-Witten curve (|5.106p on a small phase 
space. The instantonic expansion, similarly to that of the U{1) theory (|4.90p . can be developed in terms of the 
functions (|6.153p and their derivatives. For example, in [Tl we have checked, that 



•^i-Xl'^' (6.154) 



for U{2) gauge group and the only nonvanishing ti, ^2, using instantonic expansions of the equations l|5.99p . 
(jETTBl) and ([51^ . 



Elliptic uniformization for the U (2) theory 

In the case of U{1) theory the problem was solved explicitly by construction of the function (|4.7ip due to explicit 
uniformization of the rational curve (|4.69p in terms of "global" spectral parameter w. This is hardly possible for 
generic non-Abelian theory with the hyperelliptic curve (|5.99p of genus g = N — 1, but in the next to rational 
case with N — 2 

4 

2/2 = II (z- x+) = - ^0 = ^(^) (6-155) 

1=1,2 1=1 

it is an elliptic curve, and therefore can be uniformized using, for example, the Weierstrass functions 

R'{x„)/4 p'i^o) 

Z = Zq + — — — — — = Zq+ - 



^0 + C(e - eo) - C(e + Co) + 2C(eo) (6.156) 
dz ^ p'{i^)p'{0 

(p(0 - p(eo))' ^ 



where it was convenient to take 



Zq = X4 

R\zo) = a::43a;42a;4i = 4p'(^o) (6.157) 
R"{zo) = 2 {xiiXi2 + X41X43 + X42X43) = 24p(^o) 

The Abelian integral for ^{z) can be now performed in terms of the elliptic functions. Take again for simplicity 
all tk = 0, ii k > 2. Then for the differential (|5.116p one has 

2t2Z^ + siz + sn ( dz dz \ 
d<^ = — — —L—Jidz = ± 2^2^2-4 2a^ + ... (6.158) 

y z^oo \ Z Z"^ j 
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and, as was promissed before, this asymptotic fixes the coefficients 

4 

Si = -4 - ^2 y^a^z 



4 , 4 



(6.f59) 



i—l i—1 i<j 

and completely determines here the differential l|6.f 58p in terms of the curve (|6.f 55[) . For the elliptic integral 
in (|5.f f 8|) one can now write 

$ = / d^ = 2t2 (C(e + Co) + C(e - Co)) + $1 log + (6.f60) 

The constants $i_2 are easily recovered from comparison of the expansion of 

— = -2^2 (p(e - + p(e + Co)) + $1 (c(e - Co) - c(e + Co)) + $2 (6.f 6f ) 

at ^ ^ ±^0 with (|6.158p upon (|6.156p . The jumps of a multivalued Abelian integral (|6.f 60p on the elliptic curve 
(|6.f 55p . are further constrained to integers by (|5.115p and (|5.f f 9p . which can be now rewritten in the form of 
transcendental constraints for the parameters of the Weierstrass functions. 



7 Conclusion 



We have discussed in these notes the main properties of the quasiclassical hierarchy, underlying the Seiberg- 
Witten theory, which was derived in [T] directly from the microscopic setup and instanton counting. Most of 
the progress was achieved due to existence of the "oversimplified" U{1) example, naively completely trivial 
from the point of view of the Seiberg-Witten theory. However, even in this case the partition function of the 
deformed instantonic theory becomes a nontrivial function on the large phase space, being the tau-function of 
dispersionless Toda chain, and providing a direct link to the theory of the Gromov-Witten classes. The dual 
Seiberg-Witten period ("the monopole mass") satisfies the long wave limit of the equation of motion in Toda 
chain, as a function of "W-boson mass" and the (logarithm of) the scale factor. Much less transparent non- 
Abelian quasiclassical solution is nevertheless constructed using standard machinery on higher genus Riemann 
surfaces. It is also essential, that switching on higher times deform the Seiberg-Witten curve. 

The main issue now is what is going beyond the quasiclassical limit. Could at least the "simple" U{1) 
problem be solved exactly in all orders of string coupling in more or less explicit form? It is also necessary to 
stress, that the free fermion (or boson) matrix elements up to now were considered as formal series in the higher 
times, except for ti ~ log A. Their knowledge as exact functions at least of t2 could provide us an interesting 
information about physically different "phases" of the model. Another interesting and yet unsolved problem for 
the extended theory is switching on the matter by extrapolating the higher flows to tk ^ J2aLi "^a'"' what 
corresponds hypothetically to the theory with Nf fundamental hypermultiplets with corresponding masses. We 
hope to return to these problems elsewhere. 
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